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Abstract. We show the existence of continuous periodic solutions of 
the 3D incompressible Euler equations which dissipate the total kinetic 
energy. 



1. Introduction 

In what follows denotes the 3-dimensional torus, i.e. = S-*^ x x S^. 
In this note we prove the following theorem. 

Theorem 1.1. Assume e : [0, 1] — t- M is a positive smooth function. Then 
there is a continuous vector field u : x [0, 1] — )■ and a continuous scalar 
field p -.T^ X [0, 1] — > M which solve the incompressible Euler equations 

dtv + div {v iSi v) + Vp = 

(1) 

div V = 

in the sense of distributions and such that 

e{t)= \v\^{x,t)dx VtG[0,l]. (2) 



Obviously, if we choose a strictly decreasing function e(t). Theorem 11.11 
yields continuous solutions of the incompressible Euler equations which "dis- 
sipate" the total kinetic energy ^ f \v\^{x,t) dx. This is not possible for 
solutions: in that case one can multiply the first equation in ([T]) by f to 
derive 



dt^ + div + p j j = 0' 

Integrating this last identity in x we then conclude 

d r 



,x,t)dx = 0. (3) 

at Jj'i I 

Theorem 11.11 shows therefore that this formal computation cannot be justi- 
fied for distributional solutions, even if they are continuous. The pair 
in Theorem 11.11 solves ([I]) in the following sense: 



7o JT^ 



\Lp ■ V ^ Vip : V ®v + p div (p) dxdt = (4) 
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for all G C^{T^ x (0, 1);R^) and 

[ [ v-V^dxdt = Q for all V G C~(T3 X (0, 1)). 

Remark 1. In the usual definition of weak solution, is replaced by the 
same condition for divergence free test fields: therefore p disappears from the 
identity. With this alternative definition, for every weak solution v which 
belongs to a corresponding pressure field can then be recovered using 

— Ap = div div {v v) . (5) 

p is then determined up to an arbitrary function oft: this arbitrariness can 
be overcome by imposing, for instance, fp{x,t)dx = 0. However, as it is 
well-known, the equation ([5]) and the continuity of v does not guarantee the 
continuity of p. 

1.1. Onsager's Conjecture. The possibility that weak solutions might 
dissipate the total kinetic energy has been considered for a rather long time 
in the fluid dynamics literature: this phenomenon goes under the name of 
"anomalous dissipation". In fact, to our knowledge, the existence of dis- 
sipative solutions was considered for the first time by Lars Onsager in his 
famous 1949 note about statistical hydrodynamics, see [18]. In that paper 
Onsager conjectured that 

(a) C*^'" solutions are energy conservative when a > |; 

(b) There exist dissipative solutions with C*^'" regularity for any a < 5 

(note that, though Onsager's definition of "weak solution" is, strictly speak- 
ing, different from the one given above, it can be easily shown that the two 
concepts are equivalent). 

The first part of the conjecture, i.e. assertion (a), has been shown by 
Eyink in [11] and by Constantin, E and Titi in [2]. The proof of the last 
paper amounts to give a rigorous justification of the formal computation 
sketched above and leading to this is done via a suitable regularization 
of the equation and some careful commutator estimates. The second part 
of the conjecture, i.e. statement (b), is still widely open. A first result in 
that direction was the groundbreaking work of Scheffer [20] which proved 
the existence of a compactly supported nontrivial weak solution in X M. 
A different construction of the existence of a compactly supported nontrivial 
weak solution in x M was then given by Shnirelman in [21]. In both cases 
the solutions are only square summable (as a function of both space and 
time variables): it is therefore not clear whether there are intervals of time 
in which the total kinetic energy is a monotone function (indeed it is not 
even clear whether these solutions belong to the energy space L^(L^)). The 
first proof of the existence of a solution for which the total kinetic energy 
is a monotone decreasing function has been given by Shnirelman in [22j. 
Shnirelman's example is only in the energy space L°°([0, oo[, L^(M^)). 
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1.2. /i-principle. Our work [6l showed the existence of dissipative solu- 
tions for which both pressure and velocity are bounded. Besides the obvi- 
ous improvement (and the discovery of quite severe counterexamples to the 
uniqueness of admissible solutions, both for incompressible and compress- 
ible Euler), in this work we introduced a new point of view in the subject, 
highlighting connections to other counterintuitive solutions of (mainly geo- 
metric) systems of partial differential equations: in geometry these solutions 
are, according to Gromov, instances of the /i-principle. In particular the On- 
sager's Conjecture bears striking similarities with the rigidity and flexibility 
properties of isometric embeddings of Riemannian manifolds, pioneered by 
the celebrated work of Nash [T7]. Indeed, results of the same flavor as state- 
ments (a) and (b) can be proved in the case of isometric embeddings (see 
for instance [1] and the references therein): in comparing the Onsager's con- 
jecture and these results, the reader should take into account that, in this 
analogy, the velocity field of the Euler equations corresponds to the differen- 
tial of the embedding in the isometric embedding problem. All these aspects 
(and further developments for some PDEs in fluid dynamics inspired by our 
work) are surveyed in the note [8]. See also [H El [231 [251 [261 [27] . 

The understanding of Nash's construction was in a way a starting point 
for our approach to Theorem 11.11 As in the case of Nash, the solution 
of ([1]) is generated by an iteration scheme: at each stage of this iteration 
we produce an "almost solution" which solves Euler with an additional er- 
ror term. We name the resulting "perturbed" system of equations Euler- 
Reynolds system, since the error term has the typical form of the so-called 
Reynolds stress in the fluid dynamics literature (see [8] for an informal dis- 
cussion of this point). This error term converges to 0, while the sequence 
of almost solutions converge to an exact solution, uniformly in C^. At each 
stage the new approximate solution is generated from the previous one by 
adding some special perturbations, which oscillate quite fast. In our case 
the building blocks of these perturbations are Beltrami flows, a special class 
of exact oscillatory solutions to the Euler equations. Hence, the final result 
of the iteration scheme is the superposition of infinitely many (perturbed) 
and weakly interacting Beltrami flows. Curiously, the idea that turbulent 
flows can be understood as a superposition of Beltrami flows has been pro- 
posed almost 30 years ago in the fluid dynamics literature: see the work of 
Constantin and Majda [3]. Indeed, it was Peter Constantin who suggested 
to us to try Beltrami flows in a convex integration scheme. 

1.3. Comments on the proof. Though the proof of Theorem 11.11 shares 
several similarities with Nash's scheme, there are many points where our 
method departs dramatically from Nash's, due to some issues which are 
typical of the Euler equations and are not present for the isometric embed- 
dings. 
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1) Perhaps the most important new aspect of our scheme is a "transport 
term" which arises, roughly speaking, as the hnearization of the first equa- 
tion in ([T]): this term is typical of an evolution equation, whereas, instead, 
the equations for isometric embeddings are "static". At a first glance this 
transport term makes it impossible to use a scheme like the one of Nash to 
prove Theorem II. 1[ To overcome this obstruction we need to introduce a 
phase-function that acts as a kind of discrete Galilean transformation of the 
(stationary) Beltrami flows, and to introduce an "intermediate" scale along 
each iteration step on which this transformation acts. 

2) The scheme introduced by Nash and its wide generalization by Gromov 
known as convex integration heavily relies on one-dimensional oscillations 
- the simple reason being that these can be "integrated", hence the name 
convex integration. As already mentioned, the main building blocks of our 
iteration scheme are Beltrami flows, which are truly three-dimensional os- 
cillations. The issue of going beyond one-dimensional oscillations has been 
raised by Gromov (p219 of [13]) as well as Kirchheim-Miiller-Sverak (p52 
of [IS])) but as far as we know, there have been no such examples in the 
literature so far. In fact, it seems that with one-dimensional oscillations 
alone one cannot overcome the obstruction in 1). 

3) A third, more technical, new aspect is the absence of a simple po- 
tential to generate solutions of the Euler-Reynolds system: in a sense we 
cannot simply "integrate" Beltrami flows. In order to overcome this issue 
we introduce a "corrector term" to the main perturbation. This corrector 
term is not "explicit" : it is determined by solving some appropriate elliptic 
equations. 

4) In order to estimate the corrector and show that its contribution is 
negligible compared to that of the main perturbation, we use a combination 
of standard Schauder theory and oscillatory integrals estimates. This gives 
to our proof a "hard" PDE flavor compared to the construction of Nash, 
which is more on the side of "soft analysis" . 

5) As a minor comment we remark that obviously the smoothness of e in 
Theorem 1 1 . II can be relaxed, but we do not pursue this issue here. Moreover, 
the same theorem can be proved if we replace [0, 1] with [0, oo[: in this case 
we require in addition that e and its derivatives are uniformly bounded and 
that there is a positive constant cq with e > cq. 

6) In Theorem II . II our aim was to construct continuous weak solutions. 
In particular we did not address issues concerning the initial-value prob- 
lem. However, this is a typical way to proceed with problems involving 
the /i-principle. In general one may distinguish two aspects: a local and a 
global one. In geometric situations the local one is typically a differential 
constraint whereas the global one is topological (cf. jlOj). The flexibility 
(in other words the lack of uniqueness) that one observes in instances of the 
/i- principle is tied to the specifics of the local aspect. Thus, our Theorem ll.il 
deals exclusively with the "local" aspect for the Euler equations, whereas a 
possible analogue of the global aspect would be the imposition of an initial 
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data, possibly together with an admissibihty condition (as in [7]). In sub- 
sequent papers we plan to address such "global" issues (e.g. initial data, 
compactly supported ancient solutions, etc). 

1.4. Acknowledgements. We wish to thank Peter Constantin and Sergio 
Conti for several very valuable discussions on earlier attempts to prove Theo- 
rem [LTJ Moreover we are grateful to Antoine Choffrut for several comments 
on earlier versions of the paper, which considerably improved its readabil- 
ity. The first author acknowledges the support of the SFB Grant TR71, 
the second author acknowledges the support of the ERG Grant Agreement 
No. 277993 and the support of the Hausdorff Center for Mathematics in 
Bonn. 

2. Setup and plan of the paper 

The proof of Theorem 11.11 will be achieved through an iteration proce- 
dure. Along the iteration the maps will be "almost solutions" of the Euler 
equations. To measure "how far" a solenoidal field is from being a solu- 
tion of incompressible Euler we introduce a system of differential equations 
which we call Euler- Reynolds system. The name is justified by the fact that 
the matrix-field i? is a well known object in the theory of turbulence, called 
"Reynolds stress" (cp. with [8] and the references therein). In what follows 
S^^^ denotes the vector space of symmetric trace- free 3x3 matrices. 



Definition 2.1. Assume v,p,R are smooth functions on T'^x [0, 1] taking val- 

'o 



ues, respectively, m M'^, M, 5q^'^. We say that they solve the Euler- Reynolds 



system if 

dtv + div {v ®v) + Vp = div R 

(6) 

div V = . 

We are now ready to state the main proposition of this paper, of which 
Theorem 11.11 is a simple corollary. 

Proposition 2.2. Let e be as in Theorem Then there are positive 

constants rj and M with the following property. 

Let 6 < 1 be any positive number and {v,p,R) a solution of the Euler- 
Reynolds system ^ such that 

f.„)<e«)-/M^(..0 we 10, 11 P) 

and 

sup|^(x,t)| < r/<5. (8) 

x,t 

Then there is a second triple {vi,pi, Ri) which solves as well the Euler- 
Reynolds system and satisfies the following estimates: 

^e{t) < e{t) - I \vi\^{x,t)dx < ^e{t) € [0, 1] , (9) 
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sup \ Ri{x, t)\ < ^r]5 , (10) 

x,t 

sup \vi{x, t) - v{x,t)\ < mV6 (11) 

x,t 

and 

sup \pi{x, t) - p{x,t)\ < M6 . (12) 

x,t 

As already mentioned, Theorem 11.11 follows immediately from Proposition 
[221 

Proof of Theorem We start by setting vq = 0, po = 0, Rq = and 
6 := 1. We then apply Proposition 12.21 iteratively to reach a sequence 
{vn,Pm Rn) which solves ([6]) and such that 

l^<e{t)- J \vn\\x,t)dx 

sup \Rn{x,t)\ 

x,t 

SUp\Vn+l{x,t) - Vn{x,t)\ 
x,t 

sup \pn+l{x,t) -Pn{x,t)\ 

x,t 

Then {vn} and {pn} are both Cauchy sequences in C(T^ x [0, 1]) and converge 
uniformly to two continuous functions v and p. Similarly converges 
uniformly to 0. Moreover, by ()13p 

/ \v\'^{x,t)dx = e{t) ViG[0, 1]. 

Passing into the limit in ([6]) we therefore conclude that {v,p) solves ([T|). □ 

2.1. Construction of vi. The rest of the paper will be dedicated to prove 
Proposition 12.21 The construction of the map vi consists of adding two 
perturbations to v: 

Vi = V + Wo + Wc =■ V + w. (17) 

To specify the form of the perturbation Wq, which is a highly oscillatory 
function and for which we give a rather explicit formula, we need several 
ingredients. The vectorfield v + Wo is not in general divergence free. There- 
fore we add the correction Wc to restore this condition. Having added the 
correction, the main focus will then be on finding maps Ri and pi with the 
desired estimate and such that 



< foralHe[0,l] (13) 



< 



V 



2n 

< M 
M 



< 



2n 



(14) 

(15) 
(16) 



dtvi + divx{vi (8) vi) + Vpi = div^^i^i 
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The perturbation Wq will depend on two parameters, n and A, which will 
satisfy the following conditions 

A,^,-GN. (18) 

In order to achieve the estimates, A and /.i will be chosen quite large, depend- 
ing on appropriate norms of v. As already mentioned, the building blocks 
for the perturbation Wq are Beltrami flows. In order to give the formula 
leading to the definition of Wq we must, therefore, study closer the partic- 
ular "geometry" of these flows. This will be done in the next section. We 
will then be ready to define the perturbations Wo and Wc- this task will be 
accomplished in Section [5] where we also prescribe the constants t] and M 
of the estimates in Proposition 12.21 After recalling some classical Schauder 
theory in Section [5l in the Sections [6] and [7] we will prove the relevant es- 
timates of the various terms involved in the construction, in terms of the 
parameters A and fi. The choice of these parameters will be finally specified 
in Section [8l where we conclude the proof of Proposition 12.21 



3. Geometric preliminaries 

In this paper we denote by R"^", as usual, the space of n x n matrices, 
whereas 5"^" and 5q ^" denote, respectively, the corresponding subspaces of 
symmetric matrices and of trace- free symmetric matrices. The 3x3 identity 
matrix will be denoted with Id. For definitiveness we will use the matrix 
operator norm := max|^|=i \Rv\. Since we will deal with symmetric 
matrices, we have the identity \R\ = max|^|=i \Rv ■ v\. 

3.1. Beltrami flows. We start by recalling a celebrated example of sta- 
tionary periodic solutions to the 3D Euler equations, the so called Beltrami 
flows. One important fact which will play a central role in our paper is that 
the space of Beltrami flows contains linear spaces of fairly large dimension. 

Proposition 3.1 (Beltrami flows). Let Aq > 1 and let G M'^ be such that 

Ak ■k = 0, \Ak\ = ^, A_k = Ak 

for k ^1? with \k\ = Aq- Furthermore, let 

Bk = Ak + i^^x Ak£ C\ 
\k\ 

For any choice of ak & C with 11^ = a_fc the vectorfield 

W{i) = ^^kBke'"'^ (19) 

\k\=Xo 

is divergence-free and satisfies 

\W\^ 

dW{W(S)W) = V^- — (20) 
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Furthermore 

{W®W) = i W®Wdi = \ V |afc|2 fid- . (21) 
Aa 2^^^^ V 1^1 \k\) 

In other words VF(^) defined by ()19p is a stationary solution of ([1]) with 

pressure p = ——2- ^^"^ ^^^^ ^^^i^ paper we will treat the vectors Ak G 
M^, G C'^ as fixed (the choice of as prescribed in the Proposition is not 
unique, but this is immaterial for our purposes). The proof of Proposition 
13.11 is a classic in the fluid dynamics literature, but we include it for the 
reader's convenience. 



Proof. First of all observe that a^kB^k = cLkBk- Thus the vector field 
defined in ()19p is real valued. Next notice that 

divl^(0 = Y ik- Bkake'^-^ = , 
|fc|=Ao 

because k ■ Bk = for every k. 
Observe also that 

curlVF(e) = Y ikx BkUke'''^ . 
\k\=Xo 

On the other hand 

. f k , k ( k 

= Ao X Ak + Ak^ = XoBk . 

We therefore infer curl W = XqW. Since W is divergence free, div 
(W ■'V)W and we can use the well known vector identity 

\W\^ 

div {W (^W) = {W ■ V)W = VJ — L -Wx (curl W) . 



Since we have just seen that curl W and W are parallel, (|20|) follows easily. 
Finally, we compute 

W(^W = Y^ akajBk ® Sje*('=+^>« = ^ OkOjBk B^e'^^-'J^-^ . 

k,j k,j 

Averaging this identity in ^ we infer 

{W(^W)= ^ \ak\^Bk<E)Bk. 

\k\=\o 
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However, since = B_k, we get 



|fc|=Ao 



|fc|=Ao 

On the other hand, observe that the triple \^Ak,^/2-^^ x A^,-^^ forms an 

orthonormal base ofM^. Thus, 

f k , \ / k \ k k 
2Ak Afc + 2 ( — X AfcJ X Afc I + — ^ — = Id . 



This shows (j2ip and hence completes the proof. □ 

3.2. The geometric lemma. One key point of our construction is that 
the abundance of Beltrami flows allows to find several such flows v with the 
property that 

(v (g) v)(t) := ,^ ^„ / v®v{x,t)dx 

equals a prescribed symmetric matrix R. Indeed we will need to select these 
flows so as to depend smoothly on the matrix R, at least when R belongs 
to a neighborhood of the identity matrix. In view of ()2ip . such selection is 
made possible by the following Lemma. 

Lemma 3.2 (Geometric Lemma). For every iV G N we can choose rg > 
and Ao > 1 with the following property. There exist pairwise disjoint subsets 

Aj C {k e : \k\ = Xo} J G {!,..., iV} 

and smooth positive functions 

jI'^ G (Br, (Id)) i G {1, . . . , iV}, A: G A, 

such that 

(a) k G Aj implies —k G Aj and 7^"'^ = J^j,; 

(b) For each R G B^.^ (Id) we have the identity 

^ = ^ E (^I^'^ (^)) ' - ^ ^ ^) ^ (Id) . (22) 

Remark 2. Though it will not be used in the sequel, the cardinality of each 
set Aj constructed in the proof of the Lemma is indeed bounded a priori 
independently of all the other parameters. A close inspection of the proof 
shows that it gives sets with cardinality at most 98. This seems however far 
from optimal: one should be able to find sets Aj with cardinality I4. 

The proof of the Geometric Lemma is based on the following well-known 
fact. 



Proposition 3.3. The set n is dense in S^. 



10 



CAMILLO DE LELLIS AND LASZLO SZEKELYHIDI JR. 



Proof. Let s : — >• be the inverse of the stereographic projection: 
/ 2v 2u u'^ + v^ - V 

S{U,V) .- ^^^2 _^ ^2 _^ 1' ^2 _^ ^2 _^ 1' ^2 _^ ^2 _^ 

It is obvious that s(Q^) C Q"'. Since is dense in and s is a diffeomor- 
phism onto \ (0, 0, 1), the proposition fohows triviahy. □ 

Indeed, much more can be proved: ^Z^nS^, distributes uniformly on the 
sphere for re S N large whenever n = 1,2,3,4,5,6 (mod 8). This problem 
was raised by Linnik (see [IS])) who proved a first result in its direction, and 
solved thanks to a breakthrough of Iwaniec [H] in the theory of modular 
forms of half- integral weight (see, for instance, [9] and |19j). 

Proof of Lemma VJ.'A For each vector u € \ {0}, we denote by My the 
3x3 symmetric matrix given by 

V V 
= Id - -- (g) —r . 

\v\ \v\ 



R=\y. (I'ihRifMu. (23) 



With this notation the identity (j22|) reads as 

2 ^ 

fcGAj 

Step 1 Fix a Ao > 1 and for each set F <Z {k ^ "L^ : \k\ = Aq} we 
consider the set c{F) which is the interior of the convex hull, in 5^^^, of 
{Mfc : k G i^}. We claim in this step that it suffices to find a Aq and N 
disjoint subsets Fj C {fc E : = Aq} such that 

(d) -F,=F,- 

(e) c{Fj) contains a positive multiple of the identity. 

Indeed, we will show below that, if Fj satisfies (d) and (e), then we can find a 

ro > 0, a subset Tj C Fj and positive smooth functions A^:'^ G C°°(S2ro(Id)) 
such that 

(i) 

We then find Aj and the functions 7^ by 

• defining Aj := Tj U — Fj; 

• setting A^^^ = if A; G Aj \ Tj; 

• defining 



7? -74'' 

for every k £ Aj. 

Observe that the functions and the sets satisfy both (a) and (b). Moreover, 
since at least one of the A^], is positive on i?2ro(Id), 7^"'^ is smooth in i?r(,(Id). 

We now come to the existence of the set Fj . For simplicity we drop the 
subscripts. The open set c{F) contains an element aid with a > 0. Then 
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there are seven matrices Ai, . . . , Aj in c{F) such that aid belongs to the 
interior of their convex hull, which is an open convex simplex S. We choose 
•d so that the unit ball U of center aid and radius is contained in S. Then 
each point R £ U can be written in a unique way as a convex combination 
of the elements Ai : 

7 

R = ^(3,{R)Ai 

i=l 

and the functions /3j are positive and smooth on U. 

By Caratheodory's Theorem, each Ai is the convex combination ^ Xi^nMy. 
of at most 7 My. ^ with Vi^n £ F, where we require that each Xi^n is positive 
(observe that Caratheodory's Theorem guarantees the existence of 7 points 
M^. ^ such that Ai belongs to the closed convex hull of them; if we insist 
on the property that the corresponding coefficients are all positive, then we 
might be obliged to choose a number smaller than 7). 

Set ro := Then, 

R = y] -f3i{aR)Xi,nMy^ „ Vi? G S2ro(Id) . 

i,n 

and each coefficient 

-/3i{aR)Xin 
a 

is positive for every R £ i?2r()(Id). 

The set Tj is then given by {uj,n}- Note that we might have Vi^n = ^^«,m 
for two distinct pairs {i,n) and (l,m). Therefore, for k G Tj, the function 
Afc will be defined as 

Xk{R)= V -^i{aR)Xi,n. 

(i,n):k=Vi^„ 

Step 2 By Step 1, in order to prove the lemma, it suffices to find a number 
Ao and disjoint families Fi, . . . , Fjy C AqS^ H such that the sets c(Fj) 
contain all a positive multiple of the identity. By Proposition 13.31 there is 
a sequence A^ f c« such that the sets H j^"^^ converge, in the Hausdorff 
sense, to the entire sphere S^. 

Given this sequence {Afc}, we can easily partition each AfcS^ n into N 
disjoint symmetric families {-Fj^}j=i,...,Ar in such a way that, for each fixed 
j, the corresponding sequence of sets {j^F'j^}k converges in the Hausdorff 
sense to Hence, any point of c(S^) is contained in c{-^F^) provided k 
is large enough. On the other hand it is easy to see that c(§^) contains a 
multiple of the identity aid (for instance one can adapt the argument of 
Lemma 4.2 in [6]). By Step 1, this concludes the proof. □ 
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4. The maps vi, Ri and pi 



We have now all the tools to define the maps vi , Ri and pi of Proposition 
I2.2[ Recalling (fT7j) and the discussion therein, w := vi — v is the sum of two 
maps, Wo and Wc- Wq is a highly oscillatory function based on "patching 
Beltrami flows" and it will be defined first, in Section I4.1[ will then 
be added so as to ensure that vi is divergence free: in order to achieve 
this we will use the classical Leray projector, see Section Bl3] for the precise 
definition, pi is related to Wq by a simple formula, given in Section 14.41 
Finally, in Section 14.51 we will define Ri . Essentially, this last matrix field 
can also be thought of as a "corrector term", analogous to Wc- In fact, if we 
consider the point of view of [6], the Euler- Reynolds system can be stated 
equivalently as the fact that the 4x4 matrix 



is a divergence-free in space-time. Ri has therefore the same fiavor as Wc 
and is also defined through a suitable (elliptic) operator, cp. with Definition 



4.1. The perturbation Wq- We start by defining a partition of unity on 
the space of velocities, i.e. the state space. Choose two constants ci and C2 
such that < ci < C2 < 1 and we fix a function G C^{Bc2{0)) which 
is nonnegative and identically 1 on the ball i?cj(0). We next consider the 
lattice C and its quotient by {2Zf, i.e. we define the equivalence 
relation 



each /c € denote by (pk the function 

ipk{x) := ip{x - k) . 

Observe that, if k ^ i € Ci, then |A: — £| > 2 > 2c2. Hence (pk and ipi have 
disjoint supports. On the other hand, the function 




vi vi + pild — Ri vi 
vi 



1421 




is smooth, bounded and bounded away from zero. We then define 




and 
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Since a; and have disjoint supports for / / Z G Cj, it follows that for all 

l4'^K^)l' = E"'(H', (24) 

and in particular X]j=i \4'k\''^ ^ '^)\'^ — 1- Furthermore, for the same reason 
there exist for any m = 0, 1, 2, . . . constants C = C{m) such that 

suv\D'^(I)^i^\v,t)\ < C{m)fi"'. (25) 

Fix next any {v,t) and j. Observe that there is at most one I G Cj with 
the property that ai{fiv) / and this / has the property that \fj,v — l\ < 1. 
Thus, in a neighborhood of {v, r) we will have 

dr^l^^ + iik ■ v)^^^ = ik-{v-^-^ <i)^\ (26) 

Combining (j25|) and (|26|) . for any m = 0, 1,2, . . . we find constants C = 
C{m, \ k\) such that 

suplD'^idrCkl^^ + i{k-v)(l)'i^\<C{m,\k\)fi"'-\ (27) 



V 



We apply Lemma [312] with = 8 to obtain Aq > 1, > and pairwise 
disjoint families Aj together with corresponding functions ^y^^^ £ C°° {Brg (Id)). 

Next, set 

Pit) ■■=^{eit) [i-iyi^jv\\.,t)d. 

and 

R{x,t) := p{t)ld - R{x,t) , 

and define 

Wo{x,t) := yW)Yl E ivix,t),Xt)B,e^'>^-. (28) 

4.2. The constants t] and M. Note that Wq is well-defined only if — € 
i?r.p(Id) where rg is given in Lemma |3.2[ This is ensured by an appropriate 
choice of ry. Indeed, 

p{t) > ^ N3 7e(^) > niin e(i) =: c5m , 
3(27r)'^ 4 te[o,i] 

where c is a dimensional (positive) constant and m > by assumption. 
Then 
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Thus, it suffices to choose 

Observe that this choice is independent 6 > 0. 

Notice next that, by our choice of p{t) and by d?]), p{t) < 5e{t). Thus there 
exists a constant M > 1 depending only on e (in particular independent of 
5) so that 

lko||o<^^- (30) 
This fixes the choice of the constant M in Proposition 12.21 

4.3. The correction Wc- We next define the Leray projector onto divergence- 
free vectorfields with zero average. 

Definition 4.1. Let v e C°°(T^,M3) be a smooth vector field. Let 

Qv := V(/)+ -/■ V, 

where (p £ C°°(T^) is the solution of 

A(() = div V in 

with j-js (j) = 0. Furthermore, let V = L — Q be the Leray projection onto 
divergence-free fields with zero average. 

The vector field vi is then the sum of v with the Leray projection w of 
Wo, namely 

vi{x, t) := v{x, t) + Vwo{x, t) =: v{x, t) + w{x, t) 

and hence 

Wc{x,t) := —Qwo{x,t) = w{x,t) — Wo{x,t). 

4.4. The pressure pi. We define 

I I ^ 

Pi-=P (31) 

4.5. The Reynolds stress ^i. In order to specify the choice of Ri we 
introduce a new operator. 

Definition 4.2. Let v G C°°(T^,]R^) be a smooth vector field. We then 
define TZv to be the matrix-valued periodic function 

TZv := ^ {VVu + {VVuf) + ^{Vu+ (Vuf) - i(divu)Id, 



where u e C°°(T^,M'^) is the solution of 

Au = V — j' V inT"^ 

with u = 0. 



DISSIPATIVE CONTINUOUS EULER FLOWS 
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Lemma 4.3 (7^ = div"^). For any v G C°°(T^]R3) /^^^g 

(a) Tlv{x) is a symmetric trace-free matrix for each x G T^; 

(b) divTZv = V — V. 

Proof. It is obvious by inspection that TZv is symmetric. Since Vv is divergence- 
free, we obtain for the trace 

3 3 
ti{nv) = -(2divu) - -divu = 0. 

Similarly, we have 

13 1 

div [TZv) = -A(Pn) + -(Vdivn + An) - -Vdivu. (32) 

On the other hand recall that Vu = u — V(/> — j u = u — where A(j) = 
divti. Therefore A(Vu) = An — Vdivn. Plugging this identity into (|32p . we 
obtain 

div (TZv) = An 

and since n solves An = v — j- v, (b) follows readily. □ 
Then we set 

Ri := TZ {dfVi + div {vi vi) + Vpi) . 

Note that [dtVi + div (ni vi) + Vpi] has average zero. Indeed: 

• the vector field div {vi ® ni ) + Vpi has average because it is the 
divergence of the matrix field vi<^ vi -\- pild; 

• for the same reason, the identity dtv = —div (v (8> v+pld — R) shows 
that dtv has average 0; on the other hand w = Vwo has average 
because of the definition of V; this implies that dtw has also average 
zero and thus we conclude as well that dfVi = dtV + dtW has average 
zero. 

Therefore from Lemma 14.31 it follows that Ri{x,t) is symmetric and trace- 
free and that the identity 

dfVi + div {vi (8) ni) + Vpi = div Ri 

holds. 

The rest of this note is devoted to prove that the triple {vi,pi, Ri) sat- 
isfies the estimates ([9]), (jlOp . (jlip and (I12p . This will be achieved by an 
appropriate choice of the parameters fi and A. In particular, we will show 
that the estimates hold provided /i is sufficiently large and A much larger 
than /i. 

5. ScHAUDER Estimates 

In the following m = 0, 1, 2, ... , a G (0, 1), and /3 is a multiindex. We 
introduce the usual (spatial) Holder norms as follows. First of all, the supre- 
mum norm is denoted by ||/||o := supfs |/|. We define the Holder seminorms 
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as 

[/U= max WD^fh, 

\l3\=m 

[j\m+a = max sup ■ . 

m=mx^y \x-yr 

The Holder norms are then given by 

m 

= Eifh 

j=0 

\m+a — ||/||m ~l" [/]m+o- 

Recall the following elementary inequalities: 

[f]s<C{e^~V]r + e-lf\\o) (33) 

for r > s > 0, and 

[f9]r<C{[f]r\\g\\o+\\fU9]r) (34) 

for any 1 > r > 0. 

Finally, we recall the classical Schauder estimates for the Laplace operator 
and the corresponding estimates which we can infer for the various operators 
involved in our construction. 

Proposition 5.1. For any a G (0, 1) and any m S N there exists a constant 
C{a,m) with the following properties. 7/ ^ : — t- R are the unique 
solutions of 

' A<p = f 

f <A = 
then 

< C(m,a)||F|U,„. (35) 

Moreover we have the estimates 

\Qv\\m+a < C{m,a)\\v\\m+a (36) 
\'Pv\\m+a < C{m,a)\\v\\m+a (37) 

\TZv\\m+l+a < C{m,a)\\v\\m+a (38) 
||7e(div A)\[m+a < C{m, a) PlU+a (39) 
||7^Q(div A)\\m+a < C{m, a)\\A\\^+a. (40) 

Proof. The estimates ()35p are the usual Schauder estimates, see for instance 
[121 Chapter 4]. The meticulous reader will notice that the estimates in [12] 
are stated in M"" for the potential-theoretic solution of the Laplace operator. 
The periodic case is however an easy corollary. Take for instance (j) and / and 
consider them as periodic functions defined on M.^. Consider g = fx, where 
X is a cut-off function supported in Bqt^{0) and identically 1 on 54,^(0) . Let 
(j) be the potential-theoretic solution in of Acj) = g. For (J) we can invoke 




DISSIPATIVE CONTINUOUS EULER FLOWS 
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the Schauder estimates as in |12l Chapter 4] . Moreover i;^ — is an harmonic 
function in -B4pi(0). Obviously ||i;^||l2(5^^(o)) can be easily bounded using 
A(/) = /,/(/) = and the Parseval identity. Thus, standard properties of 
harmonic functions give \\4> — </'||c""."([27r]3) < (^(^-,0)11/110. 

The estimates ([36]) , ([37|) , ([38|) and ([39]) are easy consequences of ([35]) and 
the definitions of the operators. The estimate ([40p requires a little more 
care. Let u : T3 m3 ^^le unique solution of 

with j- u = 0. Then 

< C{m, a)\\A\\m+a ■ (41) 

First of all, with the argument above, one can reduce this estimate to a corre- 
sponding one for the potential-theoretic solution of the biLaplace operator 
in M^. For this case we can then invoke general estimates for elliptic k- 
homogeneous constant coefficients operators (see for instance |24] Theorem 
1]) or use the same arguments of \12\ Chapter 4] replacing the fundamental 
solution of the Laplacian with that of the biLaplacian. Finally, (|40p follows 
from the identity 

7^Q(div A) = ^ {VVu + (VVuf) + ^ (Vn + (Vuf) - ^(divn)Id 

and the estimates (|4ip and ([37p . □ 

In what follows we will use the convention that greek subscripts of Holder 
norms denote always exponents in the open interval (0, 1). 



Proposition 5.2. Let k ^I? \ {0} and \ >1 he fixed, 
(i) For any a G C°^(T^) and m (^N we have 



T3 



a{x)e 



iXk-x 



dx 



< 



(ii) Let (px G C°°(T3) be the solution of 



A^x = fx in 



with (px = 0, where 



fx{x) := a{x)e 



iXk-x 



a{y)e'^^-y dy. 



Then for any a € (0, 1) and m £ N we have the estimate 



C 



C 



c 



where C = C{a,m). 



(42) 



(43) 
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Proof. For j = 0,1, . . . define 



k ( k V' 



k 

\W 



Direct calculation shows that 



Fj{x,Xx) = — div [Aj(3;,Ax)] + —Fj^i{x, Xx). 
A A 



In particular for any m £ N 



m— 1 



a[x)e 



iXk-x 



Fo{x, >^x) = -^div [Aj{x, Ax)] + -^Fm{x, Ax) 

j=0 



Integrating this over and using that \k\ > 1 we obtain (|42p . 
Next, using psp and (j34p we have for any j < m — 1 

P,(-,A-)||a<C(A"[a], + [a],+J 

<CA^'+-(A-™H^ + ||a||o) 

and similarly 

||i^m(-,A-)|U <C(A°H^ + [a] 

m+a ) • 

Moreover, according to standard (j35p . 

||V(/>|U ^ ^(a ^ -\\A,{.,X.)\U + ^\\Fm{-A-)\\a + f^ Jo{x,Xx)dx 
hence, using ()42p for the last term, 



as required. 



□ 



Corollary 5.3. Let /c G Z'^ \ {0} be fixed. For a smooth vectorfield a G 
C°^(T3;M3) let F{x) := a(x)e*^'=-^. Then we have 

c c c 

\\'R-{F)\\a < ^^Ikllo + -J^[a]m + —[a]m+a, 

where C = C{a,m). 

Proof. This is an immediate consequence of the definition of IZ, the Schauder 
estimate ([37|) for V and Proposition 15.21 above. □ 



DISSIPATIVE CONTINUOUS EULER FLOWS 
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6. Estimates on the corrector and the energy 

In all subsequent estimates, unless otherwise stated, C denotes a generic 
constant that can vary from line to line, and depends on e, v, R as well as 
on Ao,Q! and 5, but is independent of A and /i. Smallness of the respective 
quantities will be achieved by an appropriate choice of A, ^ in Section [8l 
Moreover, all estimates will implicitly assume (jlSh . in particular that 1 < 
/i < A. 

Our aim is to estimate the space-time sup-norm oi vi — v = w = Wo + Wc 
and Ri. Since Wc and Ri are defined in terms of the singular integral 
operators Q and 7^, which act in space, instead of obtaining directly 
estimates of the norm, we will use Schauder estimates to obtain bounds 
on spatial Holder norms. Thus, in the sequel the Holder norms will denote 
spatial norms, and are understood to be uniform in time t G [0, 1]. Moreover, 
if the Holder exponent is denoted by a greek letter, then it is a number in 
the open interval (0, 1). 

It will be convenient to write Wq as 

Wo{x,t) = W{x,t, Xt, Xx), 

where 

Wiy,s,T,C):= Yl ak{y,s,T)Bke'^< (44) 
|fc|=Ao 

g 

= a/K^E E {^TV) ^"^^ «y^')'^)B,e^'-'- (45) 
i=ifceAj ^ ^ 

We summarize the main properties of the coefficients W: 

Proposition 6.1. (i) Let ak G C°°{T^ x [0,1] x M) be given by (jH]). Then 
for any r > 

\\ak{-,S,T)\\r < Cfl'', 
\\dsak{-,S,T)\\r < Cfl''^^, 
\\drak{-,S,T)\\r < Clf, 

\\{drak + i(k ■ v)ak){-,s,T)\\r < C/'"\ 
(ii) The matrix- function W ®W can be written as 

{W0W){y,s,T,O = R{y,s)+ Uk{y,s,T)e"'<, (46) 

l<|fc|<2Ao 

where the coefficients Uk G C°°it^ x [0, 1] x M;^^^^) satisfy 

Ukk = ^{tr Uk)k (47) 
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and for any r > 

\\Ui:{;S,T)\\r < Cl/, 

\\dsU>^{-,S,T)\\r < C^/' 

\\drUk{-,S,r)\\r < 

\\{drUi: + i{k.v)Ui:){;S,T)\\r < C/.^ 



In all these estimates the constant C depends on r and e, v, R hut is inde- 
pendent of (s, t) and fi. 

Proof. The estimates for are a consequence of ()25p and (j27p . Indeed, 

since the constants in the estimates are allowed to depend on e, v, R, one 

(7) 

only needs to keep track of the number of derivatives of (pj^ with respect to 
V. Then the estimates on Ok, dgOk and drak + i{k ■ v)ak immediately follow. 
Prom the triangle inequality we can then also conclude the estimate on drOk- 
Next, consider the expansion of ^ 1— t- into a Fourier series in ^, i.e. 

{W®W){y,s,T,0 = Uo{y,s,T)+ Yl Uk{y,s,T)e''<. 

l<|fc|<2Ao 

Since each Uk is the sum of finitely many terms of the form ak'Uk", the 
estimates for Uk follow from those for a^. 

Next, since Uq is given by the average (in in order to obtain (]46p we 
need to show that 

/ W^W{y,s,T,^)d^ = R{y,s). 
To this end we calculate: 

= rY, {aiiv)f = R. 

Finally, ()47p is a direct consequence of Proposition 13.11 in particular (j20p . 

□ 

After this preparation we are ready to estimate all the terms in the per- 
turbation scheme. First of all we verify that the corrector term Wc is indeed 
much smaller than the main perturbation term Wo- 



DISSIPATIVE CONTINUOUS EULER FLOWS 

Lemma 6.2 (Estimate on the corrector). 
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\\Wc\\a < Cj^ 

Proof. We start with the observation that, since k ■ = 0, 



(48) 



.\k\=\o 



k X Bk 
-JkT' 



,iXx-k 



+ 



+ y ^ iVak{x,t, Xt) 



kxB 



k ^iXx-k 



Hence 



where 



Wc{x,t) = jQUc{x,t), 



(49) 



u, 



.{x,t) = iVafc(x, t, At) 
|fc|=Ao 



The estimate (|48p then follows from the Schauder estimate (j36p for Q 
combined with 

ll^iclU ^ C/iA". 

□ 

Next, we verify the estimate on the energy Q. 
Lemma 6.3 (Estimate on the energy). 



e(t)(l - i5) - / 



< C- 



A 



(50) 



Proof. Taking the trace of identity (j46p in Proposition 16.11 we have 
\Wiy,s,T,0\^ = tTR{y,s)+ ^ Ck{y, s,T)e'''< 

l<|fe|<2Ao 

for some coefficients G C°°(T^ x [0, 1] x M), which satisfy the estimates 

\\ck{-,S,T)\\r < Cn'. 

From part (i) of Proposition 15.21 with m = 1 we deduce 



\wo\ —trRdx 



T3 



A 



and 



V ■ Wodx 



T3 



Hence, combining with (j48p we see that 
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Recalling that 

tri? = 3p = ^(e(t)(l-i<^)-|^J.Pcix), 



we conclude (j50|l . □ 
7. Estimates on the Reynolds stress 

Rewrite 

dtvi+dw {vi (S)Vi) + Vpi = [dtWo + v ■ Vwo] + 

+ div {wo<^ Wo - ^\wo\'^ld + R) (51) 

+ [dtWc + div {Vi (g) Wc + Wc (S) Vl - Wc Wc + V (g) Wo)] 

In other words we split the Reynolds stress into the three parts on the right 
hand side. We will refer to them as the transport part, the oscillation part, 
and the error. In the following we will estimate each term separately. 

Lemma 7.1 (The transport part). 

\\n{dtWo + v-Vwo)\\a<c(^^ + ^^ (52) 
Proof. Observe that 

n{dtWo + V ■ Vwo) = X] (^^"'^ + ■ ^)"'k){x, t, Xt)Bke'^^-'' | + 

\|fc|=Ao / 

+ 7^ X {dsak + v-Vyak){x,t,\t)Bke}^^-'-' 

\\k\=\o 

For the first term Corollary 15.31 with m = 2 implies the bound 

^i ^ A1-" ^ A ' 
whereas for the second term Corollary 15.31 with m = 1 implies the bound 

+ T7— T + 



Ai-° Ai-° A ■ 
Since 1 < ;U < A, we obtain ([52]) . □ 



Lemma 7.2 (The oscillation part). 



„2 

||7^(div(w;o®u;o-i|^«o|2M+i^)||„ < C-^ (53) 
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Proof. Recall the formula ()46p from Proposition 16.11 Since p is a function 
of t only, we can write the oscillation part in ()5ip as 

div {wo <8) Wo-l{\wo\'^ - p)ld + R) 

= div {wo ® Wo - R - ^{\wo\'^ - tri?)Id) 

= div iUk-^itrUk)Id)ix,t,Xt)e'^''-' 

l<|fc|<2Ao 

= Yl dwy[Uk-UtrUk)U]ix,t,Xt)e'^''-^ 

l<|fc|<2Ao 

Corollary 15.31 with m = 1 then implies (j53p . □ 



Concerning the error, we are going to treat three terms separately, as 
follows. 



Lemma 7.3 (Estimate on the error I). 

\\n{dtWc)\\a<C 



1-a 



Proof. Recah from (jM]) that Wc = jQuc. Now 



=Ao 



\k\=Xo 

Moreover, for any Ck G C°°(T'^; M"^) we have 



Ck(x,t,Xt)e'^''-'' = 4div 
iX 



k 



Ck{x,t,Xt)(S>ji^e'^^-'' 



I / k 



■ V I Ck{x,t,Xt)e 



iXx-k 



(54) 



iX \ 

Therefore dtUc can be written as 

dtUc = div Uc + Uc, 

where 

\\Uc\\a < CfiX'', \\Uc\\a < C/i^A". 

Therefore we have TZdtWc = \{TZQd\YUc + TZQuc). From the Schauder 
estimate (|40p for the operator T^Qdiv , we conclude that 

\\T^dtw^\U < J (WnQdivUcWa + WnQucWa) 

<jmc\\a + \\Uc\\a)<C^ 

□ 
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Lemma 7.4 (Estimate on the error II). 

\\R (div {vi ® Wc + Wc® Vi - Wc® Wc)) \\a < C ^l^2a ' 

Proof. We first estimate 

\\vi ® Wc+Wc ® Vi — Wc® Wc\\a < 

< CdlwillolliyclU + IbilUlkdlo + IkdlolkdU) 

m ^ 
a" 
_^ 



< C^[Z^{\\vi\\a + \\Wc\\a) 
- ^TTT^dl'^IU + ll'^clU + \\Wo\\a) 



< c-^ (c + C-^ + CX' 

Recall that 1 < < A and hence 1 < < A". Thus we conclude 

II II ^ ^ 

\\Vi ® Wc + Wc® Vi - Wc® Wc\\a < C ^i_2a ' 

(|55]) follows from the latter inequality and the Schauder estimate (|39|) . □ 



Lemma 7.5 (Estimate on the error III). 

||7^(div(^;®w;o))|U <C^. (56) 

Proof. Since Bf^ ■ k = 0, we can write 

div {v ® Wo) = Wo • V-y + (div Wo)v 

= Yl [ak{Bk-V)v + v{Bk-Vak)]e'^''-^ 
\k\=Xo 

The claim follows from Corollary 15.31 with m = 1 . □ 

8. Conclusion: Proof of Proposition 12.21 

In this section we collect the estimates from the preceding sections. For 
simplicity we set 

It should be noted, however, that due to the requirement (jlSp we can only 
ensure /U ~ A'^ for large A. 

We claim that for an appropriate choice of a and /3, the estimates (j9]l. (110p . 
(jlip and (|12p will be satisfied for sufficiently large A. First of all recall that, 
by the choice of M we have 



(57) 
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(cp. with ()30p ) and M > 1. Therefore (jlip follows from the estimate 
1 1 U^c I a < C^X"'-^ (cp. with (jll])) if, for instance, we can prescribe 

On the other hand ()12p follows easily from ()57p and the identity pi — p = 

-l\Wo\^- 

Also, from ([50]) it follows that Q is satisfied provided 

CA°+^^i < -5 min e(t). 

~ 8 te[o,i] 

Finally, (I52p . (153p as well as the estimates on the error (154p - (156p imply that 

WRiWa < C ( A"-^ + A"+2/^-^ + A2"+^-i 



tlllo 

Therefore, any choice of a, /3 such that 

Q < /3, Q + 2/3 < 1 (58) 



will ensure that (j9|l . (fTOj) . (fTT|l and (fT2]l wiU be valid for sufficiently large 
A. This completes the proof. As a side remark observe that (|58p requires 



a < i. 
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